hep-ph/ 9608^19 



CALT-68-2069 



Collision Induced Decays of Electroweak Solitons: 
Fermion Number Violation with Two Initial Particles [| 

Krishna RajagopalQ 

California Institute of Technology 
Pasadena, CA 91125, USA 

Abstract 

This paper presents work done in collaboration with E. Farhi, J. Goldstone, and 
A. Lue which is described in full in Ref. 

We consider a variant of the standard electroweak theory in which the Higgs sector 
has been modified so that there is a classically stable weak scale soliton. We explore 
fermion number violating processes which involve soliton decay. A soliton can decay 
by tunnelling under the sphaleron barrier, or the decay can be collision induced if the 
energy is sufficient for the barrier to be traversed. We discuss classical solutions to 
the Minkowski space equations of motion in which a soliton is kicked over the barrier 
by an incoming pulse. We then consider a limit in which we can reliably estimate 
the amplitude for soliton decay induced by collision with a single l^-boson. This 
amplitude depends on g like exp(— c^f"^/^), and is larger than that for spontaneous 
decay via tunnelling in the same limit. Finally we show that in soliton decays, light 
SU (2) L doublet fermions are anomalously produced. Thus we have a calculation of a 
two body process with energy above the sphaleron barrier in which fermion number 
is violated. 

1 Introduction 



In the standard electroweak theory, fermion number violation is present at the 
quantum level but these processes are seen only outside of ordinary perturbation the- 
ory. A baryon number three nucleus can decay into three leptons. The process is 
described as an instanton mediated tunnelling event leading to an amplitude which 
is suppressed by exp(— Svr^/^f^), with g ~ 0.65 the SU (2) gauge coupling constant. At 
energies above the sphaleron barrier 0], fermion number violating processes involv- 
ing two particles in the initial state are generally believed to be also exponentially 
suppressed 0]. (At energies comparable to but below the sphaleron barrier, Euclidean 
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methods have been used to show that the exponential suppression is less acute than 
at lower energies, but the approximations used fail at energies of order the barrier 
height and above.) Unsuppressed fermion number violating processes are generally 
believed to have of order 47r/(yf^ particles in both the initial and final states. This all 
suggests that fermion number violation will remain unobservable at future accelera- 
tors no matter how high the energy, whereas in the high temperature environment of 
the early universe such processes did play a significant rolef^. 

In this paper, we explore the robustness of these ideas by studying a variant of the 
standard model in which the amplitudes for certain fermion number violating colli- 
sions, as well as for spontaneous decays, can be reliably estimated for small coupling 
g. The model is the standard electroweak theory with the Higgs mass taken to infinity 
and with a Skyrme term[|^ added to the Higgs sector. With these modifications, the 
Higgs sector supports a classically stable soliton which can be interpreted as a parti- 
cle whose mass is of order the weak scale [H]. Quantum mechanically, the soliton can 



decay via barrier penetration |10|, 0. Classically, i.e., evolving in Minkowski space 
using the Euler-Lagrange equations, the soliton can be kicked over the barrier if it is 
hit with an appropriate gauge field pulse. Correspondingly, the soliton can be induced 
to decay quantum mechanically if it absorbs the right gauge field quanta. Regardless 
of whether the decay is spontaneous or induced, ordinary baryon and lepton number 
are violated in the decay. We shall see that the model has a limit in which fermion 
number violating amplitudes can be reliably estimated both for processes which occur 
by tunnelling and for those which occur in two particle collisions between a soliton 
and a single IV-boson with energy above the barrier. 

1.1 The Model 



To modify the standard model so that it supports solitons, proceed as follows. 
Note that in the absence of gauge couplings the Higgs sector can be written as a 
linear sigma model 



d^^^d^"^] - J (Tr [$^$ 



— V 



(2) 



where 

((^0; V^i) is the Higgs doublet, and v = 246 GeV. One advantage of writing the 
Lagrangian in this form is that it makes the SU {2)l x SU (2)/? invariance of the Higgs 
sector manifest. The scalar field $ can also be written as 



^ = aU 



(3) 



where U is SU{2) valued and a is a real field. In terms of these variables 



a 



(4) 



The Higgs boson mass is v2Af . We work in the limit where the Higgs mass is set to 
infinity and a is frozen at its vacuum expectation value v/^/2. Now 



= -Tr 



(5) 



which is the nonlinear sigma model with scale factor v. We will consider only those 
configurations for which the fields approach their vacuum values as |x| —>■ cxd for all 
t. We can then impose the boundary condition 



lim 

|x|^oo 



f/(x,t) 



(6) 



which means that at any fixed time U is a map from 5*^ into SU{2). These maps 
are characterized by an integer valued winding number which is conserved as the 
U field evolves continuously. However if we take a localized winding number one 
configuration and let it evolve according to the classical equations of motion obtained 
from it will shrink to zero size. To prevent this we follow Skyrme0 and add a four 
derivative term to the Lagrangian. The Skyrme term is the unique Lorentz invariant, 
SU{2)l X SU{2)ji invariant term which leads to only second order time derivatives 
in the equations of motion and contributes positively to the energy. 



^Hks ~ ^Tr 



32e2 



:Tr 



(7) 



where e is a dimensionless constant. 

Of course this Lagrangian is just a scaled up version of the Skyrme Lagrangian 
which has been used[^ |13| to treat baryons as stable solitons in the nonlinear 
sigma model theory of pions. To obtain the original Skyrme Lagrangian replace v in 
(0) t>y fn = 93 MeV. The stable soliton of this theory, the skyrmion, has a mass of 
73 /jr/e and has a size ~ 2/(e/^)[|l3l. Best fits to a variety of hadron properties give 
e = 5.5|0. The soliton of (|^) has mass 73v/e and size ~ 2/{ev) and we take e as 
a free parameter since the particles corresponding to this soliton have not yet been 
discovered. 

The standard electroweak Higgs plus gauge boson sector is obtained by gauging 
the SU{2)l X U{1)y subgroup of SU{2)l x SU{2)r in the Lagrangian (|ip. Throughout 
this paper we neglect the U{1) interactions. The complete Lagrangian we consider is 
obtained upon gauging the SU{2)l symmetry of (^: 



1 v' 
r = TrF^'' F + — 



Tr 
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Tr 



where 



D^U = id,-tA^)U (9) 

with Afj_ = 74°cr"/2 where the a" are the Pauh matrices. In the unitary gauge, U = 1, 
and the Lagrangian is 

^ = ^{-1 Tr F'^'^F^, + TiA^A^ + ^ Tr [A^ , aA , (10) 



where we have introduced 



gv , ^ 4e^ 
— and c = — 
2 



m = — and ^ = . (11) 



Note that the equations of motion derived from (|T0]) agree with those obtained by 
varying and then setting U = 1. Also note that for fixed m and ^ the classical 
equations of motion are independent of g. Since m is dimensionful and sets the 
scale, characteristics of the classical theory depend only on the single dimensionless 
parameter ^. 



1.2 The Soliton and the Sphaleron 



The classical lowest energy configuration of (|3) has = and the quantum 
theory built upon this configuration has three spin-one bosons of equal mass m. In 
the limit where g goes to zero with e and v fixed (hence, ^ goes to infinity) the 
Lagrangian is well approximated by its ungauged version (|^) which supports a 
stable soliton, so one suspects that for large ^ the Lagrangian (H) and its gauge-fixed 
equivalent (|TUp also support a soliton. In fact, Ambjorn and Rubakov[|ll| showed 
that for ^ larger than ^* = 10.35 the Lagrangian (^) does support a classically stable 
soliton whereas for ^ < ^* such a configuration is unstable. Let ?7i(x) be the winding 
number one soliton, the skyrmion, associated with the ungauged Lagrangian (|^. 
For large ^, this configuration is a good approximation to the soliton of the gauged 
Lagrangian (|]), so in the unitary gauge the soliton is ^4^°^ ~ iUldiUi, A^^^ = 0. For 
all ^ > ^* the quantum version of the theory described by (|TOD has, in addition to the 
three equal mass l^-bosons, a tower of particles which arise as quantum excitations 
about the soliton, just as the proton, neutron and delta can be viewed as quantum 
excitations about the original skyrmion |T2|, |T^. 

The Lagrangian (0) determines a potential energy functional which depends on 
the configuration A^(x). The absolute minimum of the energy functional is at A^ = 0. 
For ^ > ^* there is a local minimum at the soliton = A^°^{x.) with nonzero energy 
given by the soliton mass Mgoi- (Of course a translation or rotation of A^^^{x) produces 



a configuration with the same energy so we imagine identifying these configurations 
so that the soliton can be viewed as a single point in configuration space.) Consider a 
path in configuration space from = to A^°^ (x) . The energy functional along this 
path has a maximum which is greater than the soliton mass. As we vary the path, 
the maximum varies, and the minimum over all paths of this maximum is a static 
unstable solution to the classical equations of motion which we call the sphaleron of 
this theory. (The sphaleron of the standard modelf^ marks the lowest point on the 
barrier separating vacua with different winding numbers. Here, the sphaleron barrier 
separates the vacuum from a soliton with nonzero energy.) For fixed v and e the 
sphaleron mass Mgph goes to infinity as g goes to zero reflecting the fact that for 
g = 0, configurations of different winding (f/'s with different winding in (^) cannot 
be continuously deformed into each other. For fixed g and m, as ^ approaches ^* from 
above the sphaleron mass comes down until at ^ = ^* the sphaleron and soliton have 
equal mass. For < ^* the local minimum at nonzero energy has disappeared. 
For ^ > the classically stable soliton can decay by barrier penetration p|, [10], [11 



This process has been studied in detail by Rubakov, Stern and Tinyakov||14|] who 
computed the action of the Euclidean space solution associated with the tunnelling. 
They show that in the semi-classical limit as ^ — oo the action approaches Sir'^/g'^ 
whereas as ^ — > ^* with g fixed the action goes to zero since the barrier disappears. 

1.3 Over the Barrier 



In this paper, we focus on processes where there is enough energy to go over the 
barrier. In the standard model, the sphaleron mass is of order m/g'^ and the sphaleron 
size is of order 1/m. This means that for small g two incident W bosons each with 
energy half the sphaleron mass have wavelengths much shorter than the sphaleron 
size. This mismatch is the reason that over the barrier processes are generally believed 
to be exponentially suppressed in ly — 1^ collisions. In contrast, in the model we 
consider we can take a soliton as one of the initial state particles. To the extent that 
the soliton is close to the sphaleron we have a head start in going over the barrier. 
We can also choose parameters such that an incident W boson with enough energy 
to kick the soliton over the barrier has a wavelength comparable to both the soliton 
and sphaleron sizes. 

In Ref. |]l|, we first look in some detail at solutions to the Minkowski space 
classical equations of motion derived from the Lagrangian (0). Here, we give only a 
brief discussion of these solutions and their implications. To simplify the calculations 
we work in the spatial spherical ansatz[|l^. We solve the equations of motion for 
in the unitary gauge in the spherical ansatz numerically. We first show how to find 
the electroweak soliton for any value of ^ > Then, as initial data we take a single 
electroweak soliton at rest with a spherical pulse of gauge field, localized at a radius 
much greater than the soliton size, moving inward toward the soliton. In Ref. |jl|, we 



display one example of a soliton-destroying pulse in detail. At early times we have 
a soliton and an incident pulse and at late times we have outgoing waves only, the 
soliton having been destroyed. For ^ within about a factor of two of for all the 
pulse profiles we have tried with the pulse width comparable to the soliton size, there 
is a threshold pulse energy above which the soliton is destroyed. The energy threshold 
is larger than the barrier height, and does depend on the pulse profile.^ However, the 
existence of a threshold energy above which the soliton is destroyed seems robust, 
and in this sense the choice of a particular pulse profile is not important. 

A classical wave narrowly peaked at frequency uj with total energy E can be viewed 
as containing E/Tiuj particles. Making a mode decomposition, we can then estimate 
the number of gauge field quanta, that is W bosons, in a pulse which destroys the 
electroweak soliton. From (|10|) we see that such a pulse has an energy proportional 
to 1/(7^ for fixed m and ^. Thus, the particle number N of any such pulse goes like 
some constant over g'^. For example, at ^ = 12 the soliton-destroying pulse presented 
in detail in Ref. has g'^N ~ 2.5. At this value of ^, by varying the pulse shape we 
could reduce g'^N somewhat but we doubt that we could make it arbitrarily small. 
Upon reducing ^ towards ^* and thus lowering the energy barrier AE, smaller values of 
g'^N become possible.^ For example, at ^ = 11 we have found pulses with g'^N ~ 1. In 
the standard model, finding gauge boson pulses which traverse the sphaleron barrier 
and which have small g'^N appears to be much more challenging ||16|. Note from the 
form of ( p!OD that taking g to zero with m and ^ fixed is the semi-classical limit. In 
this limit, the soliton mass, the sphaleron mass and their difference AE all grow as 
1/(7^. The number of particles in any classical pulse which destroys the soliton also 
grows as l/g"^. 

The lesson we learn from studying classical solutions is that in the model we are 
treating, it is straightforward to find soliton destroying, sphaleron crossing, fermion 
number violating classical solutions. Particular pulse profiles are not required — 
pulses of any shape we have tried (with sizes comparable to the soliton size) destroy 
the soliton if their energy is above some shape- dependent threshold. 

"'^We have certainly not found the lowest energy or lowest particle number pulses which destroy the 
soliton. Indeed, a soliton destroying pulse with energy just above Ai? could be obtained by starting 
with a slightly perturbed sphaleron, watching it decay to the soliton, and then time reversing. We 
will see in the next section that for ^ near ^* a "pulse" so obtained would be a very long train of 
small amplitude waves, rather than a simple pulse of the kind we have used to destroy solitons in 
our numerical experiments. 

^We have worked at values of ^ ranging from 10.5 to 100. There is no reason to go to larger 
^, because as ^ becomes very large the soliton size (~ l/(m-\/J)) becomes much smaller than the 
sphaleron size 1/m) and the barrier height /S.E grows like Msoiv^- At large ^, therefore, the 
energy of soliton destroying pulses must become large compared to Msoi and large compared to the 
inverse sphaleron size. It is nevertheless a logical possibility that such pulses could be found with 
high frequencies and small values of g^N . For ^ = 50 and above, however, we have only found 
soliton destroying pulses which have large g^N . This suggests that because at large ^ the soliton is 
no longer similar to the sphaleron, we lose the advantage that we have in this model, relative to the 
standard model, in finding sphaleron crossing solutions with g^N small. 



The existence of soliton destroying classical pulses has quantum implications be- 
yond a naive estimate of the number of particles associated with a classical wave. In 
Ref. [0] we give a full and self-contained account of the relationship between classical 
solutions and the quantum tree approximation in a simple scalar theory. In a theory 
with an absolutely stable soliton, the Hilbert space of the quantized theory separates 
into sectors with a fixed number of solitons and states in different sectors have zero 
overlap |T7[] . The one soliton sector, for example, is a Fock space of states with one 
soliton and any number of mesons. The mesons are the quantized fluctuations about 
the soliton configuration and the states in the one soliton sector are scattering states 
of mesons in the presence of a soliton. No process, not even one involving large num- 
bers of mesons, connects states in the one soliton sector with states in the vacuum 
sector. In our theory, the electroweak soliton is not absolutely stable. The Hilbert 
space has sectors with a fixed number of solitons and any number of ly-bosons. How- 
ever, we argue in Ref. that the existence of classical solutions in which solitons 
are destroyed demonstrates that there are states in the zero and one soliton sectors 
of the quantum theory whose overlap in the semi-classical limit is not exponentially 
small. These states are coherent states with mean number of VT-bosons of order 1/ g'^. 
Knowing that some quantum processes exist which connect the zero and one soliton 
sectors suggests that we go beyond the semi-classical limit and look for such processes 
involving only a single incident VT-boson. 

There is an interesting limit in which we can reliably estimate amplitudes for 
single particle induced decays. Recall that for m and g fixed, as ^ approaches ^* from 
above the sphaleron mass approaches the soliton mass. We can hold m fixed and pick 
^ to be a function of g chosen so that as g goes to zero ^ approaches ^* in such a 
way that AE = Mgph — Mgoi remains fixed. We call this the fixed AE limit. It is 
different from the semi-classical limit in that as g goes to zero the classical theory 
is changing. We will argue in the next section that for ^ near ^* it is possible to 
isolate a mode of oscillation about the soliton whose frequency is near zero, which 
is in the direction of the sphaleron. This normalizable mode, which we call the A- 
mode, is coupled to a continuum of modes with frequencies uo > m. If the A-mode is 
sufficiently excited by energy transferred from the continuum modes, then the soliton 
will decay. We are able to estimate the amplitude for a single I^-boson of energy E to 
excite the A-mode enough to induce the decay. At threshold the cross section goes like 
exp(— c/(7^/^) where c is a dimensionless constant. In the same limit we can calculate 
the rate for the soliton to decay by tunnelling and we get exp(— (9/(9 — 2\/?,))c/g^^^). 
Both are exponentially small as g goes to zero and the ratio of the tunnelling rate to 
the induced decay rate is exponentially small. 



1.4 Fermion Production 



Having described classical and quantum processes in which electroweak solitons 



are destroyed, in Ref. [Q] we argue that if we couple a quantized chiral fermion to 
the gauge and Higgs fields considered in this paper, then soliton destruction implies 
nonconservation of fermion number. The argument we present treats the gauge and 
Higgs fields as classical backgrounds. In particular, we ask how many fermions are 
produced in a background given by a solution to the classical equations of motion in 
which a soliton is destroyed. We expect that our conclusions will also be valid for 
soliton destruction induced by a single VT-boson. 

We introduce fermions into this theory as in the standard electroweak theory but 
neglecting the U{1) interaction. The left-handed components transform as SU{2)l 
doublets whereas the right-handed components are singlets. The fermion mass is 
generated in a gauge-invariant way by a Yukawa coupling to the Higgs field. For 
simplicity we only consider the case where both the up and down components of 
the fermion doublet have equal mass rrif. of have the same mass rrif. The gauge 
invariant normal ordered fermion current is not conserved, that is, 

d,J^ = ^^e>^'^-^TT{F,^F^^) . (12) 

We consider backgrounds given by solutions of the kind found numerically in which 
an incoming classical pulse destroys the soliton. After the soliton has been destroyed 
the solution dissipates. By dissipation we mean that at late times the energy density 
approaches zero uniformly throughout space. This means that at late times the 
solutions are well described by solutions to the linearized equations of motion 

I A\in 



d^d" + m') = (13) 



in unitary gauge. It is tempting to try to integrate (|T^) and relate the fermion number 
violation to the topological charge 

Q = ^Jd'x e^'^^P Tr {F.^F^p) . (14) 
First, note that because the region of space-time in which F^^ 7^ is not bounded. 



there is no reason to expect Q to be an integer. Furthermore, it is shown in Ref. [0 
that for a background which satisfies ([13|) at early and/or late times the integral in 
([T^) is not absolutely convergent and Q cannot sensibly be defined. 

In a background given by a solution to the equations of motion which dissipates 
at early and late times, the number of fermions produced is known to be given by 



the change in Higgs winding number |ll8|, |T9|. In this paper, the Higgs mass is infinite 



so the Higgs winding number can never change. For solutions with no solitons in the 
initial or final states, the arguments of Ref. [0 apply, and no fermions are produced. 



However, if there is a soliton in the initial or final state the assumption of Ref. [T^ 
that the solution dissipates is not satisfied. In Ref. |T| , we show that in a background 
given by a solution in which one soliton is destroyed, one net anti-fermion is produced 
if the fermion is light (m/L <C 1 where L is the size of the soliton) and no fermions 



are produced if the fermion is heavy (m/L ^ 1). In the rrifL ^ 1 case, however, the 
sohton carries heavy fermion number [pO|, p!0| , |2T| , |l| and when the sohton is destroyed 



this quantum number is violated. In both cases there is a change of fermion number 
of minus one and heavy minus light fermion number is conserved as it must be since 
the heavy and light fermion number currents have the same anomalous divergence. 

Suppose we are only interested in light fermion production. We can view the 
heavy fermion as a device introduced only for the purpose of making an argument. 
Because we have not included the back reaction of the fermions, heavy or light, on 
the bosonic background, any conclusions we reach about the light fermion are in fact 
independent of whether there is or is not a heavy fermion in the theory. Therefore, 
in any process in which a soliton is destroyed, one net anti-fermion from each light 
SU (2) L doublet is anomalously produced. 

1.5 Relating the Model to the Real World 



The metastable electroweak soliton of the modified Higgs sector is an intriguing 
object to study. Yet this beast is not found in the standard electroweak theory 
where the Higgs sector is a linear sigma model with no higher derivative terms. It 
is reasonable to ask if the modified theory gives a credible description of physics at 
the weak scale. To date the Higgs boson has not been found. If it is found and the 
mass is low so that A of (|lD is small then working in the infinite A limit would not 
well approximate reality. However, if the Higgs is heavy, then working with infinite 
A could be justifiable. Working at the scale v and below, we then integrate out the 
heavy Higgs, leaving a low energy effective action. In this strongly interacting case, 
higher derivative terms in the effective action would not be perturbatively small and 
we would expect all possible higher derivative terms consistent with the symmetries. 
This effective theory would or would not support stable solitons. If it did then our 
use of the Skyrme term is justified as a simple way to write an effective action which 
supports solitons. 

It is possible that the Higgs is not fundamental. Rather the Higgs sector may be 
an effective theory describing the massless degrees of freedom which result 
of spontaneous symmetry breaking in some more fundamental theory. For example, 
this is the basis of technicolor theories in which the symmetry breaking is introduced 
via a scaled up version of QCD. In technicolor theories one finds technibaryons which 
can be described as electroweak solitons just as the baryons of QCD can be described 
as skyrmions. For now, regardless of whether the underlying theory is specifically a 
technicolor model, as long as we are consistent with symmetry considerations, we are 
free to choose the effective theory to conveniently describe the particles which interest 
us. Thus (0) is a simple way to describe three massless bosons (which are eaten in the 
gauged version (^) as well as a stable (metastable in (|])) heavy particle. Of course 
the effective theory includes higher derivative terms other than the Skyrme term, so 



it is not the precise form of (§) which we think is plausible, but rather the physical 
picture which it describes. 

It is worth asking what processes can sensibly be described using the effective 
theory. The effective theory is a derivative expansion in momenta over v. Consider 
the (fermion number conserving) production of soliton - anti-soliton pairs in — 
collisions. These processes are beyond the regime of applicability of the effective 
theory because the incident particles have momenta which are greater than f , and 
the underlying theory must therefore be invoked. (For example, in a technicolor 
theory the production process would be described as techniquark - anti-techniquark 
pair production followed by technihadronization.) The effective theory is, however, 
well-suited to describing soliton decay induced by a single W boson with energy just 
above AE in the fixed AE limit. In this limit m is held fixed while (7—^0, and thus 
f — > 00. Therefore, the ratio of the incident W momentum to the scale v is going to 
zero, and a treatment using the effective theory is justified. 

Over the course of this extended introduction, we have sketched all the results 
presented in full in Ref. [|I|. In the remainder of this paper, we give a complete pre- 
sentation of our treatment of quantum processes in which a single VT-boson incident 
upon the soliton kicks it over the barrier causing it to decay. In Section 3, we do a 
controlled calculation of this process in a limit in which ^ goes to ^* as g goes to zero. 
In order to do this calculation, however, we first need a better understanding of the 
classical dynamics of the theory with near and to this we now turn. 



2 Classical Dynamics for f near f 



In order to discuss the special features of the dynamics of our system for ^ near 
and because we will need it to discuss the quantum version of this theory, we 
introduce the Hamiltonian which arises from (|T^): 
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Dili' = diU'-i 
Now has no conjugate momentum and the A^ equation is 



m^A' + 



4e 



A\ A' 



A. 



+ g^D.U' = . 



(15) 



(16) 



;i7) 



This linear equation for can be solved giving A^ in terms of A^ and 11* but we do 
not need to do this explictly. The Hamiltonian for our system is given by ([T5| ) with 
A^ determined by (O) and has the general form 



H = ^IIM-\A)U+ ^V[A] 



where the sum over the coordinate x, the spatial index i and the group index are all 
implicit. The matrix M^^{A) involves derivatives with respect to x, depends on the 
configuration A, and we assume that M~^{A) is positive. Note that static solutions 
to the equations of motion, that is those with II = A = 0, occur where 6V/6A = 
and have 11 = 0. The classical equation of motion for A which arises from (|1^) is 
independent of g. Thus for the discussion of classical dynamics which we are having 
in this section, we can set g = 1. We will restore the g dependence in the next section. 

The potential energy functional V[A] has its overall scale set by m but the to- 
pography of fixed energy contours is set by ^. Ambjorn and Rubakov |TI| showed 



that for > ^* = 10.35 there is a local minimum, the soliton, whereas for < 
this minimum is absent. For ^ > ^* there is also a sphaleron, that is a saddle point 
configuration whose energy is greater than that of the soliton. As C, approaches ^* 
from above, the sphaleron and soliton merge. 

We are particularly interested in configurations which, at least initially, are small 
perturbations around the soliton. To work with these configurations we find it conve- 
nient to make a canonical transformation which has the effect of setting M~^{Asoi) = 1 



and 



dA/' 



= 0. To see that this is possible let be some complete set of orthonor- 

^sol 



dA 

mal, spatial vector, matrix-valued functions of x, indexed by a, which can be used to 
expand 11 and A. Let the coefficients of the expansion of A relative to the soliton be 
and the coefficients of the expansion of 11 be p^, that is 

/l(x,t)-A-'(x) = ^g"(t)/„(x) 

a 

n(x,t) = ^p,(t)/„(x) . (19) 

a 

(Note that the transformation from A(x, t), n(x, t) to q°'{t), Pa{t) is canonical.) Upon 
making this transformation, ([T8|) has the form 

H=^g''^{q)PaP(^ + V{q) . (20) 

A canonical transformation of the form 

q'" = q'"{q) and P'a = (21) 



can be viewed as a general coordinate transformation with pa transforming as a 
covariant vector. It is always possible to choose coordinates such that 



(22) 



is equal to 5"^ with dg'^^ / dq 
at g° 

means that the Hamiltonian 



at any given point. In fact this can be accomplished 

This 

0) can be written as 



(the soliton) with a transformation of the form g'" = C^q^ + C'^^q^q^ 
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+ 0{q') pp + V{q) 



(23) 



where we have made the required canonical transformation and dropped the primes. 
Note that V{q = 0) = M^oi and {dV / dq'^)\q=Q = 0. 

For ^ > ^* consider small oscillations about the soliton. The frequencies squared 
are given by the eigenvalues of the fluctuation matrix d'^V/dq°'dq^ at g = 0. The 
soliton is a localized object so fluctuations far from the soliton propagate freely. 
Therefore the fluctuation matrix at the soliton has a continuous spectrum above m^. 
A given soliton configuration and a translation or rotation of that configuration have 
the same energy and both solve dV/dq"' = 0. This implies that at g = there 
are six zero eigenvalues of d^V/dq°'dq^ . The associated modes which correspond to 
translating and rotating the soliton are not of interest to us and will be systematically 
ignored. 

For ^ close to ^* we now argue that there is one normalizable mode whose frequency 
ujo goes to zero as ^ goes to To see this we write 
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At the soliton (g = 0) and at the sphaleron the first derivatives are zero. As ^ 
approaches ^* the sphaleron and soliton merge so go^s to zero. It is useful to 
introduce the normalized function gsph 
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For C, > ^* the fluctuation matrix d^V/dq'^dq^ at the sohton has only positive eigen- 
values (except for the translation and rotation zero modes which play no role in this 
discussion). Equation (^) tells us that at ^ = ^* where (5 = 0, the fluctuation matrix 
has a zero eigenvalue with eigenvector gsph whereas for ^ close to ^* there is a small 
eigenvalue, tUg, whose associated eigenvector is close to gsph- Note that gsph points 
from the soliton to the sphaleron. Thus the low frequency mode, which we call the 
A-mode, is an oscillation about the soliton close to the direction of the sphaleron. 

For ^ > at the sphaleron there is one negative mode, that is one negative 
eigenvalue of the appropriately defined fluctuation matrix. As ^ comes down to ^* 
the sphaleron and soliton become the same configuration so this negative eigenvalue 
must come up to zero in order for the spectra of the fluctuation matrices of the soliton 
and sphaleron to agree at ^ = ^* . Therefore for ^ close to ^* the unstable direction 
off the sphaleron has a small negative curvature. There are two directions down from 
the sphaleron. One heads toward the soliton and the other heads (ultimately) to the 
classical vacuum at A = 0. We see that for near ^* the soliton can be destroyed by 
imparting enough energy to the A-mode since it is this mode which is pointed towards 
the sphaleron and beyond. 

We wish to describe the interaction of the A-mode with the other degrees of 
freedom. We use the Hamiltonian written in the form (^3]). At this point it is 
convenient to make an orthogonal transformation on the {q'"} so that the transformed 
set are the eigenvectors of the soliton fluctuation matrix d^V / dq°'dq^\q=Q. We will 
label these vectors as q^ where u"^ is the eigenvalue of the fluctuation matrix. The 
eigenfunctions include: 

i) The continuum states with eigenvalues u'^ > rn?. (Note that for each 
cu^, in general, there is more than one eigenvector. The extra labels on q^ 
are suppressed in our compact notation.) 

ii) The normalizable state q^^^ = A with eigenvalue oj^ which goes to zero as 
^ goes to i*. 

iii) The zero eigenvalue states associated with translation and rotation. 

iv) Other normalizable states which might exist but whose frequencies do not 
have any reason to approach zero as ^ goes to ^* . 



Up to cubic order the Hamiltonian ( ^31) is 
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where in the elhpses we now include all terms with modes of type iii) and iv) as 
well as higher order interactions of the A-mode and the continuum modes, p is the 
momentum conjugate to A and is the momentum conjugate to g^. The number b 
and the functions c, d and e are determined by the soliton configuration. For example 
d{uj) is presumably peaked at values of u which correspond to wavelengths of order 
the size of the soliton. As ^ goes to ^* we know that Uq goes to zero but we expect 
no dramatic behavior of b, c, d ot e in this limit. 
Consider the A-mode potential 

V{X) = \ulX' + h' + .... (29) 

There is a local minimum at A = which is the soliton and a local maximum at 
A = —uj^/b, which is approximately the sphaleron, where the second derivative is 
— cJq. We work with ^ sufficiently close to ^* that uq is small. This means that A 
at the sphaleron is small and if we only study dynamics up to and just beyond the 
sphaleron we are justified in neglecting the quartic and higher terms in A. We also see 
that as ^ goes to ^* so that Uq goes to zero, the soliton and sphaleron come together 
and at ^ = ^* the A potential has an infiection point at A = and the soliton is no 
longer classically stable. 

In order to discover the relationship between uq and (^ — ^*) as ^ approaches 
it is necessary to study the behavior of the A-mode potential as ^ approaches In 
( ppp for every value of ^, we have shifted A so that the minimum of the potential is at 
A = 0. This ^ dependent change of variables obscures the behavior of the coefficients 
of the potential before the shift. Calling the unshifted variable A, then if we expand 
the potential in terms of e = ^ — ^* about e = where there is an inflection point, we 
have 

V{\, e) = 0(e)A + 0{e)\^ + (b + 0(e)) A=^ + . . . , (30) 

where 6 is a constant. We know that the coefficients of A and A^ are zero at e = 0, 
and we assume that these coefficients can be expanded about e = and we know of 
no reason for the order e terms to vanish. For e > the minimum of the potential is 
at A ~ e^/^, (A is shifted relative to A by this amount), and at the minimum of the 
potential d^V/d\^ ~ e^/^ that is 

A small amplitude oscillation of the A mode will decay because of its coupling to 
the continuum modes which can carry energy away from the soliton. However for 
uq < m this decay is very slow in the sense that the characteristic time for the decay 
is much greater than 1/ujq. To understand this consider A(t) as a source for radiation 
in the continuum via the coupling J^dujd{u!)q^ in the Hamiltonian (^). Suppose 
that X{t) is a purely sinusoidal oscillation with frequency and with an amplitude 
which is small. Radiation with frequency uq is not possible because the continuum 



frequencies begin a.t u = m. However, has frequency 2ujo and therefore if cuq > rn/2 
the couphng will excite propagating modes with uj = 2ujq and the A oscillation will 
radiate at twice its fundamental frequency. Because the coupling is of order A^, the 
rate of energy loss will be small. If cuq < m/2 then radiation at a; = 2ujq is also 
not possible. However, ifm/3 < ujq < m/2 the X^q^j coupling (which we have not 
written in ( pSD because it is fourth order) allows the A oscillation to radiate at three 
times its fundamental frequency. There is another source of radiation with u = 3ujq. 
The potential for the A-mode is not exactly quadratic so the A oscillation, although 
periodic, is not exactly sinusoidal. If the period of the oscillation is 2tt/ujo, A will be 
a sum of terms of the form sinwot, sm2u}ot, sinSuot,... with diminishing coefficients. 
This means that A^ will also be a sum of terms of this form. Those terms in A^ 
with frequencies greater than m will excite radiation via the A^g^j coupling. As ujq is 
reduced from m toward zero, the radiation is produced only by higher order couplings 
and by higher harmonics, and therefore the amplitude is reduced and the decay takes 
longer. 

We have numerical evidence for this behavior within the spherical ansatz. In 
unitary gauge in the spherical ansatz, is written in terms of four real functions of 
r and t. Gauss' law specifies one of these functions in terms of the other three, so that 
a solution to the equations of motion is fully specified by the real functions p{r,t), 
9{r, t) and ai(r, t) defined in Ref. To watch an oscillating soliton radiate for a long 
time, we implement energy absorbing boundary conditions at the large r boundary 
of the simulation lattice, as described in Ref. fl]]. We wish to excite the A-mode and 
watch it oscillate. We describe in Ref. a convenient way of choosing a configuration 
which is a soliton plus a small perturbation where the perturbation is preferentially in 
modes with lower frequencies. We then use this configuration as the initial condition 
for the equations of motion. The resulting evolution is shown in Figures 1 and 2 for 
^ = 10.4. The functions p, 6, and ai (we show p only) all oscillate about the value 
they take at the soliton and the period of oscillation is 16.69 m~^. We identify this 
with the A-mode and so obtain uq = 0.3764 m. Furthermore we see that away from 
the soliton there is a small amplitude train of outgoing radiation. After a brief initial 
period during which any perturbations not in the A-mode radiate away, the outgoing 
radiation settles down to a frequency 1.129 m, three times the fundamental frequency. 
(At r = 10 m~^, we see in Figure 2 that the frequency Suq oscillation of p has a small 
modulation with frequency uq. This is the tail of the A-mode oscillation and is not seen 
at larger values of r.) The radiation causes the amplitude of the A-mode to decrease 
very slowly — by about 4% over 80 oscillations. We have done similar simulations 
at ^ = 11 and ^ = 12 also, where we find ujq = 0.80 m and ujq = 0.98 m respectively. 
In these simulations, the oscillating soliton emits radiation with u = 2ujo, and the 
amplitude of the radiation and the rate of decay of the fundamental oscillation are 
larger than in Figure 2. The values of Uq for ^ = 10.4, 11, and 12 which we have 
found numerically are in good agreement with the relationship (^Tj). This numerical 
evidence suggests that we are justified in using the Hamiltonian (^) to describe the 
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Figure 1: As described in the text, we have perturbed the ^ = 10.4 sohton and let it 
evolve for a long time. Here, we show p(r) for a series of different times: t = 0, 144, 
288, . . . 1440 m^^. This shows the envelope of the oscillation of p. In Figure 2, we 
show p at r = 0.608/m and r = 10/m as a function of time. 

long-lived normalizable A-mode with cuq < m and its coupling to the continuum. In 
the next section we will quantize this Hamiltonian and use it to describe the excitation 
of the A-mode by single VF-boson quanta. 

Finally we note that in principle it is possible to destroy a soliton with a minimum 
energy pulse, i.e. one whose energy is just above AE, and for ^ close to ^* this energy 
is small. To find the form of this pulse we could time reverse a solution which starts 
at the sphaleron and is given a gentle push towards the sohton. For ^ close to ^* so 
that the A-mode has a small frequency, the configuration takes a very long time to 
settle down to the soliton and in the process emits a very long train of low amplitude 
outgoing waves. Although the time-reversed solution consisting of a very long train 
of incoming low amplitude waves being absorbed by the soliton would eventually go 
over the sphaleron barrier and result in soliton decay, it would be rather difficult 
to set up initial conditions which produce this complicated, finely tuned, incoming 
configuration. Thus, the minimum energy soliton destroying pulses are not easy to 
build although we have seen that with some extra energy, for ^ near the soliton is 
easily killed. 
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Figure 2: In the left panels, we show p at r = 0.608/m as a function of time. It 
oscillates with period 16.69/m, and the amplitude of the oscillation is decreasing very 
slowly. In the right panels, we show p at r = 10/m, to display the outgoing travelling 
waves shed by the oscillating soliton. These waves have three times the frequency 
of the fundamental oscillation seen at r = 0.608/m. Note that the amplitude of 
the outgoing waves is so small that they arc invisible in the plots of p(r) on the 
preceding page. We conclude that for ^ = 10.4 the soliton has an almost stable mode 
of oscillation with frequency ujq — 0.374 m — the A- mode — which slowly radiates 
waves with frequency Scuq. 



3 Quantum Processes in the Fixed AE Limit 



In the previous section we saw that for ^ close to ^* it is possible to identify a low 
frequency vibration of the soliton, the A-mode, with frequency much less than m. 
If enough energy is transferred to this mode the soliton will decay. In this section we 
discuss the quantum mechanics of this mode. In this quantum setting the soliton can 
decay by barrier penetration as well as by being kicked over the barrier by a single 
ly-boson. We will see that if we work in a limit where AE is held fixed as we take g 
to zero, then we can reliably estimate the leading terms in both the tunnelling and 
induced decay rates. 

The Hamiltonian for just the A-mode coming from p8| ) is given by 

where we have restored the g dependence. Note that uq, b and all the terms in the 
ellipses depend on ^ and m but not on g. We have changed the sign of A for later 
convenience. As ^ goes to ujq goes to zero but the other terms are presumed not 
to change much. The classical soliton is at A = while the sphaleron is at A = c^Q/fo 
from which we have 

^^-Iw- '''' 

The fixed AE limit has g going to zero with ^ taken to ^* in such a way that (|33|) 
is fixed. Since b{^,m) does not vary much as ^ goes to we see that in this limit 
^0 ~ 9^^^- Using (1^) and (pID, we see that g'^AE ~ — so that in order to 

take the fixed AE limit we take g to zero with — ^*) ~ g^^^. (The reader who is 
concerned that the coefficient of A^ in (|3^), uj^/g"^, goes to infinity in the fixed AE 
limit should note that because of the g^ in front of the in (^) the frequency of 
oscillation is uq.) 

When taking the fixed AE limit, it proves convenient to rescale according to 

A' = Xuo/gr^Xg-'/' 

h' = bg/oulr^bg^ (34) 

Writing the Hamiltonian ( p2D in terms of the new variables and then dropping the 
primes we obtain 

H, = ujI^ + V{X) , (35) 



where 

V{X) = 1a^ - ^A^ + . . . . (36) 



After rescaling, the sphaleron is at A = 1 /6 and the barrier height is given by 



AE = 1/(66") . (37) 

Quartic and higher terms in V^(A) are all suppressed by powers of (yf/cjo ~ g'^''^- Note 
that ujq now plays the role of h in the Hamiltonian (|35D. As g goes to zero in the 
fixed AE limit, c^o goes to zero like g^^^ and a semi-classical (WKB) treatment is 
appropriate in order to compute the leading small-^f behavior of the soliton destruction 
cross-section. 

In the fixed AE limit, the ground state of the quantum soliton has the A degree 
of freedom in a wave function ipo{^) which is described approximately by a harmonic 
oscillator ground state wave function: 

/ 1 X 1/4 / \ 2 \ 

XTTUJo/ \ 2uJo J 

There are three relevant scales in A, which differ in their (yf- dependence. First, the 



width of the ground state wave function ^ {ipolX^lipo) goes like ycJo ~ g^^^- The 
second scale, which goes like g^, is the distance in A between the sphaleron at A = 1/6 
and the minimum at A = 0. Note also that ( p8D is a good approximation to ipQ for 
A such that the cubic term in V^(A) can be neglected relative to the quadratic term, 
namely for |A| -C 1/6. Finally, note that the quartic and higher terms in V{X) can 
be neglected for A less than of order Uo/g ~ g'"^^^, the third scale. Hence, as g is 
taken to zero with AE fixed, truncating the potential at cubic order becomes valid 
for larger and larger A. 

The soliton will decay if the A degree of freedom tunnels under the barrier given 
by the potential V{X) shown in Figure 3. The rate is of the form 



Ce-"^ (39) 



where the factor B is 



B^^r'dAv^^v^r^^l^^^^. (40) 

luq Jo ^ 5 (^06 5 Uq 

We are able to neglect the width of the wave function (|38|) in this calculation because 
as g goes to zero it is small compared to the change in A during the tunnelling 
process. Since in the fixed AE limit cuq ~ g^^^ we see that the tunnelling rate goes 
as exp(— constant/fy'^/^). For the approximation to be reliable we require that B be 
much greater than one. This in turn requires that g be small. 



We can compare this calculation with that of Rubakov, Stern and Tinyakov|n 
who numerically calculated the action of the Euclidean space solution which tunnels 
under the barrier. They used the equations of motion of the full 3-1-1 dimensional 
theory with the restriction to the spherical ansatz. At ^ = 12 we have AE = 1.2 m/g"^, 




Figure 3: The potential V{X) for real A. For later use, the energies Eq and E are also 
shown. tpQ has three turning points, and A = Aq is the left-most of the three. ipE has 
one turning point at A = Ag. 



= 0.98 m giving g^B = 4.4 which is to be compared with what we read off Figure 
2 of Ref. |T^, namely g^B = 4± 1. This agreement again supports the view that the 
A mode is the relevant degree of freedom for discussing soliton decay for C, near 

We now turn to induced soliton decay. Our picture is that the soliton will decay 
if the A-mode is excited to a state with energy above AE. The A-mode couples to 
the continuum modes which can bring energy from afar to the soliton. The free 
quantum Hamiltonian for the q^j is 



where 
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(41) 



(42) 



The q^j have been chosen to diagonalize the fluctuation matrix at the soliton. There- 
fore Hq^ describes non-interacting massive IV-bosons propagating in a fixed soliton 



background. For each value of uj there are actually an infinite number of different 
W^-boson quanta. For example there are the states with frequency uj and all values 
of angular momentum relative to the soliton center. These extra labels are omitted 
throughout but their presence is understood. 

The A-mode couples to the continuum modes through cubic couplings of the form 



Hint = ~ I / '^^ '^('-^) + I duj' 



e{uj,uj')q^q^i\ (43) 



which appear in (|28|). We have rescaled A according to (^). The couplings ( ^Sf ) arose 
upon expanding about the soliton. The functions d{uj) and eiu^oj') are peaked at 
values of uj corresponding to wavelengths of order the size of the soliton. They are 
also only peaked if the unspecified labels allow large overlap with the soliton. For 
example even with oj chosen so that (tu^ — m?)^^^'^ ~ soliton size, it is only the low 
partial waves which have d{uj) and e{uj^uj') large. 

The first term in ( ^3]) allows for the absorption of a single M^-boson by the soliton. 
The VT-boson energy E is transferred to the A-mode. The second term in ( ^3] ) allows 
a single VT-boson to scatter inelastically off the soliton, transferring energy E to the 
A-mode. We now calculate the rate for the absorption process; the calculation for the 
scattering process is similar. (The coefficients of the A and A^ operators have different 
5f-dependence, but this will not affect the leading gf-dependence of the cross-section 
for either process.) Assuming that the soliton starts in its ground state, in order for 
the soliton to decay we require E + ujq/2 > AE. Since Uq -C AE we can approximate 
this as E > AE. In the fixed AE limit we are free to choose AE to be a constant 
times m where the constant is of order unity. (Recall that m is held fixed throughout 
this paper.) Now the soliton size is roughly 2/{^m^/^) and in the fixed AE limit ^ goes 
to ^* = 10.35. Thus the VT-boson wavelength and the soliton size can be comparable. 
There is no length scale mismatch and d{E) need not be small. 

Using Fermi's Golden Rule we now calculate the cross section for W + soliton — > 
anything with no soliton. Let |k) be a single W^-boson state with energy E, normalized 



to unit particle fiux. Now 

X' f „, „. A^ 
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where we have defined d(k) so that it is independent of g (see (|4^)). The A-mode 
starts in the state V'o(A) with energy ~ cc;o/2 which again we neglect relative to AE. 
The interaction (|4^) can cause a transition to a state iPe{X) in which the A-mode has 
energy E. Since the width of tpQ is ~ g^^^ ^ 1, it is tempting to try approximating 
the states with E > AE as plane waves 



The cross section for a transition from ipo to ipE is 



C^destruction = ^f (^^) ^(^)' ' (^6) 

where A/" is a (^-independent constant and where I{E) is the integral 

I{E) = I rfA^o(A)AVs(A) . (47) 

If we take iJjq and ^/^e as in (|38|) and (|^) respectively, is easily evaluated, yielding 

liE) ~ exp i-E/iOo) , (48) 

where we have dropped all prefactors. This result is in fact incorrect While it is true 
that (^) and (^5]) yield a good approximation to the integrand where the integrand 
is biggest, the result ( |48| ) is exponentially smaller than the integrand. This raises 
the possibility that corrections to the wave functions neglected to this point may 



change (^81) . We must, therefore, use WKB wave functions which take into account 
the quadratic and cubic terms in the potential V{X). As (7 in the fixed AE 
limit, ujq ^ and using semi-classical wave functions becomes a better and better 
approximation. We show below that for E = AE the leading dependence of the of 
X(£') as (7 in the fixed AE limit is in fact that of (^8[) with the coefficient of 
AE/ujq being (18 — 4-\/3)/5 instead of 1. Thus, we will find that even though the 
soliton destruction process does not involve tunnelling, the correct cross-section is 
exponentially small as c^o ~ g^^^ goes to zero. The reader who is not interested in the 
details of the evaluation of T{E) can safely skip to equation (0) . 
We now wish to evaluate the leading semi-classical dependence of 

IiEo,E) = JdXtPEo>^'^E (49) 

in the fixed AE limit where E > AE and AE > Eq > and where i/je and ipEo ai^e 
WKB wave functions for the Hamiltonian (P3|). See Figure 3. The reader may be 
concerned that (^) is infinite. (Both wave functions are real, and for large positive 
A the integrand (^) has a non-oscillatory piece which grows like A^A~^/^.) However, 
when the relevant limits are taken correctly, the answer we seek is in fact finite. 
Recall that our problem reduces to that of the A mode in a cubic potential only 
for |A| < too/g ~ g~'^^^- Therefore, we should do the A integration from A = —A 
to A = +A, where A is real and positive and where we take A to infinity more 

slowly than g'"^^^ as g goes to zero. The result of such an evaluation would go like 

^3/2 

exp (—const ant /ciJo). Because we do not take A to infinity before taking g to zero, 
the prefactor does not make the result infinite. 

■^We are grateful to D. T. Son for noticing this, and for pointing us toward the correct answer. 





Figure 4: This figure is a sketch showing the important points in the complex A-plane 
discussed in the text. The branch points are marked with dots, and the branch cuts 
are shaded. The contour in (5.19), just above the real axis, and the deformed contour 
we use to evaluate the integral are both shown. 



The evaluation of matrix elements of operators between semi-classical states has 
been treated by Landau p2|, and although his final answer does not apply to our 
problem, we follow his method to its penultimate step. Landau's method yields 
only the leading {i.e. exponential) dependence of such matrix elements, and says 
nothing about the prefactors. Thus, using Landau's method yields the leading small- 
g dependence of (^) irrespective of whether the prefactors make the integral infinite. 
In the calculation which follows, it nevertheless proves convenient to multiply the 
integrand in (^) by exp(— JA^/cuq) with J a constant. This does in fact render 
the integral finite, but it may also modify the exponential dependence of the result. 
Therefore, after the g ^ limit has been taken we must take the J — > limit. 
Landau's method applied to our problem yields 
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In this expression, A is treated as complex and it is understood that the contour 
has been deformed into the upper half plane. This is done both in order to avoid 
the turning points on the real axis shown in Figure 4, and because in deriving (|50|) 
Landau uses expressions for WKB wave functions which are valid only in the upper 
half plane and not on the real axis. The first square root in the exponent in (|50| ) is 
taken to be positive on the real axis for A < Aq and the second is taken to be positive 
on the real axis for A < A^;. 

The equation V{x) — E = has three roots. One is at Ag, on the negative real 
axis, and the other two, at Abp and A^p, have nonzero imaginary parts. (For E — > AE, 
Abp goes to the real axis at Agph = !/&•) In evaluating ( pOD we must keep in mind that 
at A = Abp in the upper half plane, the integrand has a branch point. This singularity 
will play an important role in our analysis. (Unlike in the example treated explicitly 
by Landau, it does not arise from a singularity in ^(A).) The branch cut from Abp 
must not cross the real axis, and it is convenient to take it to run upward vertically. 
The integrand in (|50|) is a function which is analytic in the upper half plane except 
at Abp and along the associated cut. To evaluate the integral, we are free to push the 
contour upward away from the real long as we ensure that it does not touch 

the branch point Abp or cross the branch cut. 

We now evaluate the leading exponential dependence of (|50|).f] To this end, we 
drop the prefactors in (|^). We write the integral as 



f dXexp — \x{X) + iY{X)] (51) 



where X and Y are real and where 

X + tY = dx ^2 {V{x) - Eq) - r dx ^2 {y{x) - E) - . (52) 

J Xq J Xe 

It is easy to check that for J = the integrand in (^TJ) has no saddle points at finite 
A. However, making J nonzero introduces a saddle point at large |A| which moves off 
to infinity as J is taken to zero and it is convenient for us to evaluate the integral 
with nonzero J and then take the J — limit. 

We now describe the behavior of X at large |A|. Write A = Aexpi^. We have 
chosen the branch cut to run vertically and so it is at 6' = 7r/2 for large A. To the 
right of the cut, that is for 7r/2>^>0, asA goes to infinity 

X ~ -A^/2 sin(5^/2) - JA^ cos(20) , (53) 

and the J term is subleading. X goes to +oo at large A for 7r/2>^^>27r/5 and goes 
to — oo for 27r/5 > > 0. The descent to — oo is most rapid for 6* = 7r/5. To the left 
of the cut, that is for tt > > 7r/2, as A goes to infinity 

X ~ X* + A-i/2 sin(^/2) - JA^ cos{2e) , (54) 



*The analysis described below and the result ( pO| ) were provided by A. V. Matytsin. 



where X* is a constant independent of J, A, and 6. (For J = 0, as A goes to infinity 
for Tc > 9 > 71/2, X X* and Y —>■ 0.) For nonzero J, there is a saddle point at 
finite A. For small J, this saddle point is at ^ ~ 37r/5 and A ~ J~^/^. Thus, as J — 
the saddle point recedes to infinity as promised, and JA' at the saddle point goes to 
zero. Therefore, in the J ^ limit X at the saddle point goes to the value X*. 

We now deform the contour as sketched in Figure 4. For nonzero J, the saddle 
point is at finite A and we choose the contour to follow the path of steepest descent 
from this saddle point. To the left of the saddle point, the steepest descent path 
curves toward the real axis, and then approaches the real axis asymptotically. As we 
discuss below, X(Abp) is greater than X*. Therefore, to the right of the saddle point, 
the path of steepest descent from the saddle point cannot get around the branch point 
and necessarily runs into the branch cut. After reaching the cut, the next section of 
the path ascends as it traverses (II), following the cut inward toward the origin, until 
it reaches the region of the branch point Abp. Along (II), X ascends monotonically 
from X* to X(Abp). Y is not constant. Then, to the right of the cut, the contour 
follows the path of steepest descent (III) toward infinity along 6 = 7r/5. 

There are two contributions to the integral (|5T|). First, the saddle point makes a 
contribution which goes like exp(X*/co'o). (Note that we take the g ^ limit and 
then take the J — > limit.) The second contribution arises because the path must 
ascend from the saddle point at infinity as it traverses (II) in order to get around 
the branch point, before then descending along (III) to the right. Therefore, the 
integral (|5lD receives a contribution from the region of the branch point which goes 
like exp(X(Abp)/u;o). In sum, therefore, the integral ( |5D| ) goes like 

J(Eo, E) ~ exp(X7^o) + exp(X(Abp)/a;o) , (55) 

where we have dropped the prefactors, about which Landau's method says nothing. 
At this point, we can take the Eq ^ and E AE limits simply by setting 
£^0 = and E = AE. Prior to this point in the calculation, taking these limits would 
require careful treatment of branch points. Henceforth we set Eq = and compute 
I{E) =I{0,E). 

It only remains to evaluate the relative size of X(Abp) and X*. Both X(Abp) and 
X* depend on E. After some calculation one finds that for E = AE 

X* = -ujqB = -—AE , (56) 
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where B is the tunnelling amplitude computed in (|40|), and 

(57) 

SO X(Abp) is the larger {i.e. least negative) of the two at E = AE. At large E, 
both X(Abp) and X* decrease like — For E > AE, the integrals in ( p2D must 



be evaluated numerically. We find that both X(Abp) and X* decrease monotonically 
with increasing energy, and X(Abp) is always greater than X* . Consequently, the 
integral is dominated by the region of the branch point for all energies E > AE. 
That is, 

J(E) ~ exp(X(Abp)M) (58) 

and 

CTdestruction ~ (k) CXp (2X( Abp) /t^o) , (59) 

where we have dropped all prefactors except d. Thus, although the integrand has a 
saddle point (at infinity), the integral is not dominated by that saddle point. This 
occurs because the path of steepest descent from the saddle point necessarily runs 
into the branch cut. Equivalently, the presence of the branch cut prevents the actual 
contour of integration from being deformed into a path of steepest descent through 
the saddle point. Although the path can be deformed to pass through the saddle, it 
must ascend from the saddle to the region of the branch point. (Note that although 
-^(Abp) > X* for all energies E > AE, X(Abp) is greater than Buq, and the rate for 
induced soliton decay is greater than the tunnelling rate, only for E within a range 
of energies which we determine numerically to be AE < E ^ l.TAAE.) 

Because X(Abp) decreases monotonically with increasing E, the cross section (|59|) 
for the soliton to be destroyed by a single VF-boson is least suppressed by X(i?) at 
threshold. For E = AE the soliton destruction cross section goes like 



^^AE/uq] (60) 

as — > in the fixed AE limit. 

We expect d{E) and accordingly d{h) to be appreciable when E ~ AE so long as 
AE is comparable to the inverse soliton size, which is of order the inverse W^-mass. 
Under these conditions, there will be no length scale mismatch and d{E) will not 
depend sensitively on E for E ~ AE, so cTdestruction will be maximized for E = AE. 
Thus the maximum rate for soliton decay induced by collision with a single VT-boson 
is proportional to exp(— (36/5 — 8\/3/5)AE/uq). This is to be compared with the 
tunnelling rate in the same limit which is proportional to exp(— (36/5)Aii^/u;o). Both 
go to zero as g goes to zero like exp(— constant/(7^/^), but the ratio of the tunnelling 
rate to the induced decay rate is exponentially small. 

We have computed the cross section for a single py-boson to be absorbed by the 
soliton and to excite the A-mode to a continuum state above the barrier, which in 
our picture results in soliton decay. The cross section for a VT-boson to destroy the 
soliton by scattering off the soliton and transferring energy E to the A-mode can be 
calculated using the second term in (^3]). The calculation is similar to the one we 
have done and the result has the same exponential factor as in (^0]) but would have 
a different prefactor. Because the exponent in (|^) includes cUq^ ~ g~^^^, these cross 
sections go to zero faster than any power of g a.s g goes to zero in the fixed AE limit. 



Note that this suppression arises even though the process does not involve tunneUing 
and even though there is no length scale mismatch. It arises as a consequence of the 
limit in which we have done the computation, because in that limit destroying the 
soliton reduces to exciting a single degree of freedom to an energy level infinitely many 
(~ AE/uq) levels above its ground state. Thus, taking — > at fixed AE makes the 
computation tractable but makes the induced decay rate exponentially small, albeit 
larger than the tunnelling rate. 

4 Concluding Remarks 

We have described a theory which agrees with the standard electroweak model 
at presently accessible energies but which includes a metastable soliton with mass of 
order several TeV. This Higgs sector soliton may have a dual description as a bound 
state particle made of more fundamental constituents or it may be that the Higgs 
sector is fundamental and when quantum effects arc taken into account, a metastable 
soliton is found. In any event, given the soliton, under certain circumstances we can 
reliably estimate the rate for collision induced decays. The parameters of the theory 
can be chosen so that the soliton configuration is close to the sphaleron configuration, 
which means that using the soliton as an initial particle makes it easy to find sphaleron 
crossing processes. Indeed, we have found classical solutions in which the soliton is 
destroyed where the incoming pulse corresponds to a quantum coherent state with 
~ VF-bosons. The rate for such processes is not exponentially suppressed as 
g goes to zero. Furthermore in the limit g goes to zero with AE — Mgph — ^soi 
fixed we can reliably estimate the rate for a two particle scattering process in which 
a single incident W^-boson kicks the soliton over the barrier causing it to decay. We 
have argued that in all processes in which the soliton disappears fermion number is 
violated. This model may be relevant only as a theoretical foil, as a demonstration 
that fermion number violating high energy scattering processes can be very different 
than in the standard model. However if no light Higgs boson is discovered, it is even 
possible that Nature may be described by such a model. 
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